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Abstract 

Geobacter species are of great interest for environmental and biotechnology applications as they can carry out direct 
electron transfer to insoluble metals or other microorganisms and have the ability to assimilate inorganic carbon. Here, we 
report on the capability and key enabling metabolic machinery of Geobacter metallireducens GS-1 5 to carry out CO2 fixation 
and direct electron transfer to iron. An updated metabolic reconstruction was generated, growth screens on targeted 
conditions of interest were performed, and constraint-based analysis was utilized to characterize and evaluate critical 
pathways and reactions in G. metallireducens. The novel capability of G. metallireducens to grow autotrophically with 
formate and Fe(lll) was predicted and subsequently validated In vivo. Additionally, the energetic cost of transferring 
electrons to an external electron acceptor was determined through analysis of growth experiments carried out using three 
different electron acceptors (Fe(lll), nitrate, and fumarate) by systematically isolating and examining different parts of the 
electron transport chain. The updated reconstruction will serve as a knowledgebase for understanding and engineering 
Geobacter and similar species. 
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Introduction 

Microorganisms play a major role in the global carbon cycle. 
Insights into the various mechanisms and energetic constrains 
which govern their behavior wiU advance our understanding of 
carbon fluxes and might ultimately allow for a rational perturba- 
tion of the carbon cycle. Key features of the carbon cycle are the 
conversion of organic and inorganic carbon and the energy flow 
through the system. Qualitative and quantitative knowledge of 
carbon assimilation and electron flow to and from key microor- 
ganisms is critical when evaluating certain aspects of the carbon 
cycle. Members of the genus Geobacter are ubiquitous in the soil 
environment, and have been described to utilize various organic 
substrates while transferring electrons to insoluble metals exter- 
nally [1]. Furthermore, certain Geobacter species, such as G. 
metallireducens have been reported to transfer electrons directly to 
poised electrodes [2] and even to other microbes [3], a process 
coined direct interspecies electron transfer or DIET. Quantitative 
assessment of carbon and energy flow in G. metallireducens by 
computational modeling approaches therefore provides valuable 
insight into the role of this bacterium in the carbon cycle. 

Geobacter metallireducens, the first Geobacter species that was isolated 
[4] , serves as a pure culture model for the study of many of the 
important reactions that Geobacter species catalyze in the 



biogeochemistry of anaerobic soils and sediments, groundwater 
bioremediation, and several bioenergy applications [5]. For 
example, Geobacter species play a major role in the biogeochemical 
cycling due to their ability to couple the oxidation of organic 
compounds to the reduction of Fe(III) and Mn(IV) oxides [5] . G. 
metallireducens was the first microorganism shown to be capable of 
the anaerobic degradation of aromatic hydrocarbons [6] and 
Fe(III) is an important electron acceptor for the removal of 
aromatic hydrocarbons in many contaminated subsurface envi- 
ronments [5]. Direct electron transfer from electrodes to 
microorganisms to drive anaerobic respiration has potential 
appfications in bioenergy and bioremediation [7]. 

Constraint-based reconstruction and analysis (COBRA) is a 
powerful method for characterizing the content of an organism, or 
systems of organisms, and understanding the limits of its collective 
functionality [8]. Metabolic network reconstruction (the most 
widely utilized form of COBRA) enables the enumeration of the 
genome-wide machinery (i.e., enzymes, uptake systems, etc.) in an 
organized fashion for use in modeling [9]. With a reconstructed 
network for an organism, predictions can be made about its 
functionality when combined with physiological data in a 
modeling framework. Further, a validated and accurate network 
can be utilized for prospective design and engineering of cellular 
networks [10]. 
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Author Summary 

The ability of microorganisms to exchange electrons directly 
with their environment has large implications for our 
knowledge of industrial and environmental processes. For 
decades, it has been known that microbes can use 
electrodes as electron acceptors in microbial fuel cell 
settings. Geobacter metallireducens has been one of the 
model organisms for characterizing microbe-electrode 
interactions as well as environmental processes such as 
bioremediation. Here, we significantly expand the knowl- 
edge of metabolism and energetics of this model organism 
by employing constraint-based metabolic modeling. 
Through this analysis, we build the metabolic pathways 
necessary for carbon fixation, a desirable property for 
industrial chemical production. We further discover a novel 
growth condition which enables the characterization of 
autotrophic (i.e., carbon-fixing) metabolism in Geobacter. 
Importantly, our systems-level modeling approach helped 
elucidate the key metabolic pathways and the energetic 
cost associated with extracellular electron transfer. This 
model can be applied to characterize and engineer the 
metabolism and electron transfer capabilities of Geobacter 
for biotechnological applications. 

There is a history of modeling Geobacter sp. using COBRA [1 1]. 
One of the first studies utilizing COBRA and Geobacter was for G. 
sulfwreducens [12]. The key findings of this study were an initial 
reconstruction and examination of the extracellular electron 
transport, the examination of the efficiency of internal biomass 
biosynthetic pathways, and predictions of gene deletion pheno- 
types. A subsequent study branched off to build a reconstruction of 
G. metallireducens based on the original content of the G. suljurreducens 
reconstruction [13]. The initial G. metallireducens reconstruction was 
used to examine the efficiency of pathway usage in the network 
along with yield on a variety of substrates. 

In this work, an updated reconstruction was built and analyzed 
to better understand key capabilities of G. metallireducens. The 
updated reconstruction effort was fueled by the appearance of a 
more complete genome annotation [14] and new data available 
for the electron transport chain and key metabolic content [15,16]. 

Results/Discussion 

Construction, comparison, and validation of the 
Geobacter metallireducens GS-15 genome-scale metabolic 
reconstruction 

An updated reconstruction of G. metallireducens GS-15, iAF987, 
was generated by reconciling an existing genome-scale reconstruc- 
tion [1 3] and an updated genome annotation, performing a bottom- 
up reconstruction of additional metabolic pathways. This new 
reconstruction was functionally tested for performance under 
known growth conditions (Figure lA). The final reconstruction 
contained 987 genes, 1284 reactions, and 1109 metabolites. In the 
first phase, the existing reconstruction was compared to the updated 
genome annotation [14] to identify a list of agreements, discrep- 
ancies, and scope for expansion. A distinct periplasm compartment 
was determined to be important as G. metallireducens has the unique 
ability to transfer electrons extraceUularly [5] . Thus, characterizing 
the electron transfer pathways from the cytosol through the 
periplasm to the extracellular space was crucial for understanding 
this unique capability. Furthermore, the addition of the periplasm 
compartment allows for a more accurate representation of 
metabolism, such as ^-cresol and 4-hydroxybenzyl alcohol 



degradation, which partially occurs in the periplasm [17]. A wild- 
type and a reduced 'core' biomass objective function [18] were 
formulated to validate whether the reconstruction could generate 
the appropriate biomass components necessary to rephcate, and for 
use in simulation to predict the growth rates of the organism on the 
diflferent substrates. Gaps were filled in the network using data 
characterizing growth of G. metallireducens GS-15 on 19 different 
carbon sources/electron donors with Fe(III) as the electron 
acceptor, and the SMILEY algorithm [19] (see Text SI). 

The AF987 reconstruction was compared to the previous 
version [13] and an automatically generated reconstruction from 
the ModelSEED framework [20] (Figure IB) was used to identify 
and evaluate newly reconstructed and unique content. The 
ModelSEED reconstruction was found to have 1 1 4 unique genes 
that were not present in the iAF987 reconstruction. Of these, 86 
genes were involved in macromolecular synthesis, DNA replica- 
tion, and protein modffications that are beyond the scope of a 
metabolic network, and 8 of them did not have a specific reaction 
association in the ModelSEED (i.e., generic terms such as 
aminopeptidase, amidohydrolase). Of the remaining 20 genes, 
only two (Gmet_0988 and Gmet_2683) were added to the 
reconstruction as isozymes for existing reactions; the other 18 
assignments conflicted with our functional annotation of the 
genome and thus were not included. The AF987 reconstruction 
contains 227 genes not in either reconstruction, thus representing 
a significant advancement of coverage. These newly included 
genes encode several unique pathways, encoding 325 unique 
reactions, which have not previously appeared in a collection of 1 4 
representative reconstructions (Table 1) from the UCSD database 
from which iAF987 was constructed and is internally consistent 
(see Text SI for a detailed comparison of AF987 to previous 
work). 

Transcriptomic data profiling a growth shift from acetate to the 
aromatic compound benzoate was integrated with the metabolic 
model to validate its content. Specifically, the computational 
analysis was performed using the MADE algorithm [21] which 
uses the statistical significance of changes in gene expression to 
create a functional metabolic model that most accurately 
recapitulates the expression dynamics. Of the 987 genes in the 
metabolic model, transcriptomic data indicated that the expres- 
sions of 857 genes do not change significantiy and 130 were 
differentially expressed (>2-fold and p-value<0.05). Specifically, 
77 genes were up-regulated and 53 were down-regulated during 
this shift. The MADE algorithm predicted that during this 
metabolic shift, the expression of 885 genes in iAF987 do not 
change significantly and 102 genes are differentially expressed. Of 
the 102 differentially expressed genes, the MADE algorithm 
predicted the up-regulation of 70 genes and down-regulation of 32 
genes. Specifically, the model predicted upregulation of 70 genes, 
while data indicated 77 genes to be upregxilated during this shift. 
Similarly while the model predicted downregulation of 32 genes, 
the data actually indicated that 53 genes were downregulated 
during this shift. The model-based prediction of change in 
expression disagreed with the in vitro transcriptomic data for only 
28 of the 987 genes leading to 97% overall agreement (for a more 
detailed breakdown, see Text SI and Table SI). Among the genes 
differentially expressed during the shift, the genes encoding for 
benzoyl-CoA reductase were up-regulated over 100-fold during 
benzoate growth. It was determined that this key enzyme that links 
the degradation of aromatic substrates to central metabolism is not 
ATP driven as previously thought [13], but is likely membrane 
bound and proton translocating [16]. Thus, a proton translocating 
reaction was added to the reconstruction for this step in 
metabolism. A translocation stoichiometry of 3 protons per 
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A. Reconstruction Workflow 



B. Comparison of Reconstruction Content 



Existing Geobacter 
metallireducens Reconstruction 
• G: 747 • R: 697 • M: 769 



Remove Lipid, Membrane, IVIurein, 
and LPS metabolism content 



List of Agreements 



Updated Genome Annotation 
• Genes: 3630 



List of Discrepancies 



Compartmentalization 
of content: cytosol, '' 
periplasm, extracellular 



Curated 3 Compartment G. 
metallireducens Reconstruction 
• G: 974 • R: 1164 • IVl: 1096 



Manual Curatlon 



KEGG 
Database 

Publication 
& Literature 
Bibliome 




Model SEED 
I Database 

MetaCyc 
Database 



Bottom-Up Reconstruction 



Final Geobacter 
metallireducens Reconstruction 
•G:987 -R: 1284 • M: 1109 



Biomass Objective 
Function for G. metallireducens 
& Culture Growth Data for 
Known Conditions 




178 



, 227 



34 



529 



53 



108 



2009 G. metallireducens 
reconstruction (iJS747) 

Model SEED 



2013 G. metallireducens 
reconstruction (iAF987) 



C. Validation of Network Content 

Experiment 



Growth 

Substrate 1 ' 

Q Substrate Shift 

-O 



Anal 



Growth 
Substrate 2 




Analysis of Network and Data 

Expression Profile 
Each Growth 
Substrate 
Condition 



iAF987 model & 
MADE algorithm 
predicted 
\ expresson change |/ ^\| 



significantly 
differentially 
expressed genes 



Agreement of 97% between Measured Values & 
Feasible Expression Transitions in the Model 



Figure 1. The workflow developed to generate iAF987 along with comparison and validation of its content. A) The workflow detailing 
the reconstruction process of the G. metallireducens metabolic network. The reconstruction process was initiated by comparing the updated genome 
annotation for G. metallireducens to the existing reconstruction to create a list of discrepancies that was manually reviewed and curated. Content that 
was in agreement with the updated annotation and reconstruction was used to generate a draft set of intracellular reactions. Lipid, membrane, murein, 
and LPS content were removed from this list as a periplasm compartment was added to the reconstruction. Manual curatlon [9] was aided by the KEGG 
[46], ModelSEED [20], and MetaCyc [35] databases. Further, numerous publications and literature sources (i.e., the "bibliome") were used to refine the 
network content. The manual review process resulted in a draft reconstruction that was used in conjunction with a formulated biomass objective 
function in simulations to validate the content of the reconstruction and generate a final version. Some figure images adapted from [47]. B) Venn 
Diagram showing the comparative analysis of gene content included in different versions of a G. metallireducens reconstruction. C) A schematic of the 
validation of network content with transcriptomics data for a shift from acetate growth conditions to benzoate with the MADE computational algorithm. 
doi:1 0.1 371 /journal.pcbi.1 003575.g001 



Table 1. Subsystem distribution of reactions unique to the G. 
metallireducens iAF987 reconstruction. 





Subsystem 


Number of Reactions 


Transport 


75 


Lipids and Glycan metabolism 


64 


Vitamins & Cofactor Biosynthesis 


51 


Aromatic Compound degradation 


26 


Alternate Carbon Metabolism 


21 


Amino Acid Metabolism 


18 


Energy Metabolism 


17 


Central Metabolism 


15 


Nitrogen and Sulfur Metabolism 


12 


Nucleotide Metabolism 


11 


Other 


9 


Metal Respiration 


6 


Total 


325 



doi:l 0.1 371/journal.pcbi. 1003575.1001 



electron was determined to be the likely extent of coupling through 
a thermodynamic analysis (see Text SI). Similar transcriptomic 
analyses for growth shifts on two other aromatic electron donors 
(i.e., toluene and phenol) yielded 86% and 84% agreement, 
respectively (Table SI). These findings will likely broaden our 
knowledge of how G. metallireducens can be utilized for bioremedia- 
tion. 

Analysis of unique metabolic capabilities: Carbon fixation 
and external electron transfer pathways 

The genome of G. metallireducens GS-15 encodes two out of the 
six known carbon fixation pathways [22]. The pathways which 
were reconstructed in iAF987 are the reductive citric acid cycle 
(rTCA) and the dicarboxylate-hydroxybutyrate cycle [22] 
(Figure 2A). Key enzymes for the rTCA include the 2-oxoglutarate 
synthase (abbreviated OOR2r in the reconstruction) and ATP- 
citrate lyase (ACITL), both of which enable the citric acid cycle to 
run in reverse. For the dicarboxylate-hydroxybutyrate cycle, the 
key enzyme is 4-hydroxybutyryl-CoA dehydratase (4HBCOAH). 
The rTCA and the dicarboxylate-hydroxybutyrate cycles share 
four reactions. The reconstruction of these carbon frxation 
pathways led to the prediction of a new growth condition for G. 
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metallireducens (Table 2). With the expanded content including the 
carbon fixation pathways, it was computationally predicted that G. 
metallireducens can grow with formate as the electron donor and 
Fe(III) as the electron acceptor when a computational screen of all 
possible media combinations was performed with the model 
(Table 2). Investigating the resulting flux distributions revealed 
that the CO2 derived from formate oxidation is reduced via the 
rTCA to form acetyl-CoA which is subsequently assimilated into 
biomass. The electrons derived from formate oxidation are split 
between running the rTCA and for Fe(III) reduction. The energy 
gained by formation of a proton gradient during Fe(III) reduction 
was instrumental for providing the required ATP for carbon 
fixation. This prediction of growth on formate and Fe(III) was 
experimentally vaUdated, Figure 2B. The requirement of CO2 
fixation for growth of G. metallireducens solely on formate and CO2 
as carbon sources is further highlighted by a study which examined 
G. sulfurreducens reducing Fe(III) with formate as the electron donor. 
While G. sulfiirreducens was able to reduce Fe(III) with formate as 
electron donor, it required the addition of 0.1 mM acetate to 
assimilate cell carbon (i.e., grow) [23]. This was attributed to the 
lack of an rTCA in G. sulfurreducens, specifically due to the absence 
of the ATP-dependent citrate lyase. Overall, this example 
represents a unique power of the model to rapidly generate 
hypotheses in silico that can then be verified experimentally. This 
particular result proves to be of great interest for examining 
carbon fixation. 

Examining ttie cost of external electron transfer 

To reconstruct the electron transport system of G. metallireducens, 
key steps involving electron transfer to the terminal electron 
acceptor were subject to a thermodynamic analysis (Text SI). 
Specifically, the feasibility- of proton translocation and the theoret- 
ical maximum proton translocation stoichiometry was determined 
[24,25]. Subsequendy, the included content was analyzed with 
physiological data from growth screens to predict stoichiometrics for 
the key reactions of the electron transport system (ETS). This was 
performed in a model-driven iterative process with an approach that 
delineates the energetics of extracellular electron transfer by 
examining three distinct modules (Figure 3A— B). These modules 
were chararterized by representative electron acceptors; fumarate, 
nitrate, and Fe(ni), respectively. Reduction of fumarate in G. 
metallireducens rec^ms the strain to harbor the dcuB gene [26]; nitrate 
is reduced via the ETS and the nitrate reductase [14,27]; 
extracellular Fe([II) requires several extracellular c-type cyto- 
chromes and pili [28], which have been shown in the closely 
related G. sulfurreducens to have metal-like conductivity [29]. Acetate 
consumption has previously been described as a preferred optimal 
growth condition for G. metallireducens [27], therefore it was used as 
the primary electron donor when analyzing the ETS in this study. 
The pathways and key ETS reactions for this conversion are shown 
in Figure 3C and evidence for the inclusion of the reactions in 
iAF987 is given in the Text SI. 

The first step in the modehng process was to examine the 
fumarate reductase reaction and maintenance energies necessary 
for predicting phenotypes using constraint-based analysis (Step 1, 
Figure 3C). To examine this content, data was generated using the 
G. metallireducens dcuB strain and utilized in simulations. An 
electrogenic bifunctional fumarate reductase/succinate dehydro- 
genase (FRD2rpp) was included based on recent biochemical 
evidence in a similar species [15]. This inclusion was validated 
with the growth data of G. metallireducens on acetate and fumarate 
using the dcuB strain. The earlier version of the G. metallireducens 
model [13] and the current model (iAF987) without an electro- 
genic fumarate reductase were unable to produce a feasible 



solution when constrained with experimentally measured acetate 
and fumarate uptake rates (see Table 3). However, the electrogenic 
fumarate reductase enabled the model to reproduce the experi- 
mental obser\'ations and was assigned to translocate two protons 
per two electrons transferred to fumarate. It was determined that 
two protons were necessary to drive the endergonic oxidation of 
succinate with menaquinol as part of the TCA cycle during Fe(III) 
respiration and nitrate reduction (see Text SI). Maintenance 
energies were also evaluated in this step by choosing a growth 
condition that eliminated any use of the external electron 
transferring reactiouM (mo these' redactions could be examined in 
isolation later). G. metallireducens dcuB strain was grown in batch 
using acetate/fumarate medium and also in a chemostat at a set 
growth rate (Table 3). The growth rate determined in this study 
was similar to that previously reported (0.1 14 hr ' vs 0.105 hr 
respectively) [26]. From these two conditions, it was possible to 
estimate the maintenance costs (GAM and NGAM) for the model 
using an established procedure [9,18]. The calculated costs were 
79.20 mmol ATP gDW"' for the GAM and 0.81 mmol ATP 
gDW"' hr"' for the NGAM. These values are similar to those 
found previously for G. sulfurreducens [12], thus it provided 
confidence in the use of the GAM and NGAM throughout the 
study. The additional content included in the ETS was further 
analyzed using experimental data. 

The second step in the modeling process to examine extracel- 
lular transfer was to examine the energetics of the menaquinone 
cytochrome oxidoreductase reaction (CYTMQPR3) in the ETS 
(Step 2, Figure 3B). This was performed by switching the electron 
acceptor from fumarate to nitrate, while keeping acetate as the 
electron donor. Assuming the nitrate reductase translocates two 
protons per two electrons (a value that has been verified [30]), the 
translocation stoichiometry of the menaquinone cytochrome 
oxidoreductase can be isolated. Analyzing the energetics of the 
CYTMQPR3 reaction, it was determined that a likely number of 
three protons are translocated per two electrons transferred to the 
cytochrome pool (see Text SI). Using this stoichiometry and 
simulating growth with acetate and limiting nitrate (see Table 3), 
the predicted optimal growth rate was calculated to be 0.054 hr ' 
as compared to the experimentally determined value of 
0.050 hr~' (in this experiment, the growth rate is equal to the 
set dilution rate in a steady state chemostat). Thus, it was 
concluded that the proton translocating assignment of the 
CYTMQPR3 reaction was consistent with the experimental 
results of the growth screen. 

The final step taken to predict the energetic cost of transferring 
electrons to an external substrate was to examine growth of G. 
metallireducens on acetate and Fe(III) after reconciling the other 
components of the ETS consistent with observed phenotypic data. 
Figure 3C, step 3 shows the reconstructed path from the inner 
membrane cytochromes to the outer membrane cytochromes and 
eventually to reduce Fe(III). An energetic: cost was estimated in 
terms of an ATP cost proportional to the flux of electrons to Fe(III). 
Growth screens of wild type G. metallireducens (acetate/Fe(III)) were 
performed in triphcate in both batch and chemostat cultures. When 
culturing with Fe(III) as an electron acceptor, a measurement of 
biomass is challenging given that the optical density is used to 
calculate the amount of reduced iron (see Methods). Therefore, the 
ratio of acceptor produced to donor consumed was used to compare 
to simulations, as it is more precise than a biomass-normalized 
uptake and production rate that was calculated by measuring 
protein content in the chemostat. Furthermore, the ratio of acceptor 
to donor calculated for both the batch and chemostat conditions was 
very similar, thus providing a consistent experimental comparison 
(see Table 3). A Phenotypic Phase Plane (PhPP) analysis [31] was 
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A. Carbon Dioxide Fixation Pathways B. Validation of modeling prediction: 

VACOAi , 4H BC0AH A novel autotrophic growth condition 




Figure 2. Reconstruction, analysis, and wet lab validation of carbon fixation pathways in G. metallireducens. A) A map of the two carbon 
fixation pathways included in the iAF987 reconstruction and encoded by genes annotated in the updated genome annotation. The two pathways, 
the reductive citric acid (TCA), and the dicarboxylate-hydroxybutyrate cycles share four reactions in the citric acid cycle. CO2 fixation (and resulting 
acetyl-coA generated) and ATP-driven steps are shown in green and red, respectively. B) A graph of phenotypic data that demonstrated the model- 
predicted growth condition of formate as an electron donor and carbon source with Fe(lll) as an electron acceptor is a viable growth condition. 
Carbon fixation is occurring under these conditions and the specifics of this process can be further elucidated in subsequent "drill-down" studies. 
doi:10.1371/journal.pcbi.1003575.g002 



used to compare model-predicted performance to the experimen- 
tally measured acceptor to donor ratio. It was calculated that a cost 
of one proton translocated across the inner membrane per one 
electron transferred ultimately to Fe(III) best matched the line of 
optimality in the PhPP analysis (see Figure 3 C, Figure SI). Further, 
this cost is very close to 0.3 ATP per electron transferred ultimately 
to Fe(III), as the ATP synthase in the cell converts protons to ATP at 
a ratio of 3.33 protons per ATP (see Text SI). Further analysis of 
this modeling approach with phenotypic data on different electron 
donors (butanol, ethanol, and pyruvate) yielded the same cost of 
external electron transfer (see Table S3). Thus, it was hypothesized 
that this is the approximate cost for external electron transfer to iron 
and the reactions and costs were built into the iAF987 reconstruc- 
tion as such. This cost can now be further validated for diflFerent 
external electron transfer processes that G. metallireducens is known to 
carry out. 



Conclusion 

The work presented here demonstrates how constraint-based 
modeling and reconstruction can be applied to generate hypotheses 
that can be tested experimentally. Specifically, modeling revealed a 
non-obvious culturing condition where carbon fixation could be 
directly examined. Further, the cost of external electron transfer 
could be quantified using an iterative and systematic approach. 
Carbon sequestration is of great biotechnological interest [32]. By 
understanding the mode of growth for CO2 fixation, computational 
predictions can be used to guide genetic modifications which enhance 
the rate of CO2 fixation. Specifically, this could be in the form of 
reaction knockouts or identification of genes which could be targeted 
for overexpression which are precficted to enhance CO2 frxation. The 
reconstructed model advances our knowledge for this unique species 
and provides a platform for further analysis and hypothesis 
formulation for environmental and biotechnology applications. 



Table 2. Validated and predicted carbon sources and electron donors for G. metallireducens. 





Experimentally and Computationally Validated Carbon 
Sources and Electron Donors 


Computationally Predicted Carbon 
Sources and Electron Donors 


Experimentally and Computationally 
Validated Electron Acceptors 


3- Methylbutanoic acid, p-Cresol, 4-Hydroxy-benzyl alcohol, 

4- Hydroxybenzaldehyde, 4-Hydroxybenzoate, Acetate, Butanol, 
Butyrate (n-C4:0), Benzoate, Benzaldehyde, Benzyl alcohol, 
Ethanol, Isobutyrate, Phenol, Propionate (n-C3:0), Propanol, 
Pentanoate, Pyruvate, Toluene 


Formate', Glycerol, Long Chain Fatty 
Acids (Saturated and Unsaturated) 


Fe(lll), Manganese Mn(IV), Uranium U(VI), 
Technetium Tc(VII), Vanadium V(V), Nitrate, 
Nitrite 


^Formate was validated as a carbon source and electron donor experimentally in this study. 



doi:l 0.1 371 /journal.pcbi.1 003575.t002 
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A. Iterative Process Utilized for Predicting 
the Cost of Electron Transfer 
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The Electron Transport System 
in G. metallireducens 

A 3 Step Process to Examine the 
Cost of External Electron Transfer 



Figure 3. A model-driven analysis of the electron transport system in G. metallireducens. A) schematic of the iterative loop process 
common to model-driven analyses, and (B) the process applied to examine the cost of the electron transport system in G. metallireducens. C) Map of 
the updated ETS in the iAF987 reconstruction. The functional states of the network components during internal electron transfer to fumarate are 
shown in green (facilitated by the dcuB strain), additional components active during internal transfer to nitrate are shown in blue, and additional 
components during external electron transfer (Fe respiration) are shown in brown. Note that the fumarate reductase (FRD2rpp) operates in opposite 
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directions depending on whether the electron acceptor is fumarate or nitrate/Fe(lll). Also note that nitrate reductase (N03R3pp) is not used when 
Fe(lll) is the electron acceptor. Abbreviations are defined in Dataset S2. The process was started with the bottom-up reconstruction of the updated 
metabolic network. At each cycle around the loop, optimal performance was calculated as different components of the ETS were isolated, compared 
to experimental data, and then incorporated into the reconstruction. Ultimately, the final product is a reconstructed ETS consistent with experimental 
data and an estimate of the cost associated with transferring electrons from the internal membrane cytochrome (focytC) to the extracellular electron 
acceptor Fe(lll) (fe3). 
doi:1 0.1 371/journal.pcbi.1 003575.g003 



Methods 

Updated annotation 

Curation of the genome annotation of G. metallireducens was 
continued after the initial publication [14] with additional insights 
from curation of the genome annotations of Geobacter bemidjiensis 
[33], Pelobacter carbinolicus [34], and other species (M. AHujkar, 
unpublished), with extensive reference to the MetaCyc database 
[35], The updated annotation was submitted to NCBI with 
reference numbers CP000148 and CP000149. 

Reconstruction process 

The reconstruction was generated in a four-step process. First, 
the updated genome annotation for G. metallireducens was entered 
into the UCSD SimPheny (Genomatica, San Diego, CA) database 
[14]. Next, the existing G. metallireducens reconstruction [13] was 
mapped to the updated genome using the existing gene-protein- 
reaction associations [8] and all of the pathways excluding 
membrane lipid biosynthesis, lipopolysaccharide biosynthesis, 
murein biosynthesis and degradation, and transport were entered 
into the SimPheny framework if an exact match for the reaction 
was present in the UCSD SimPheny database. If an exact match 
for the reaction did not exist in the UCSD SimPheny database on 
the level of metabolites participating in the reaction, they were 
manually evaluated for inclusion (see below). Next, a comparison 
of the metabohc content included in the updated genome 
annotation (Dataset SI) that was not in SimPheny was performed. 
Manual evaluation of new content or disagreements from the 
annotation and existing genome-scale reconstruction consisted of 



gathering genetic, biochemical, sequence, and physiological data 
and reconciling this information to determine the likelihood of 
each reaction being present in the organism. This manual curation 
process has been described and reviewed several times [8,9]. In the 
manual review process, the KEGG database (www.genome.jp/ 
kegg/), the ModelSEED database [20], and primary literature (see 
Dataset S2) were used extensively in the manual curation process. 
Confidence scores were given for each reaction along with 
noteworthy evidence used to justify inclusion of a given reaction. 

Reconciling the reconstruction content with 
ModelSEED. The gene content in the G. metallireducens recon- 
struction was compared and reconciled with the automated 
reconstruction obtained from ModelSEED [20] . The ModelSEED 
reconstruction was found to have 119 unique genes that were not 
present in the G. metallireducens reconstruction. Upon further analysis, 
it was determined that 86 of these genes were involved in 
macromolecular synthesis, DNA replication, and protein modifica- 
tions that are beyond the scope of a metabolic network. Of the 
remaining 33 genes, eight did not have a specific reaction association 
in the ModelSEED (i.e., generic terms such as aminopeptidase, 
amidohydrolase). The ModelSEED annotations of the 25 metabolic 
genes were compared with the updated genome annotation 
presented in tliis work. It was found that these two sets of annotations 
were consistent for only two genes (G. metaJlireducenset_0988 and G. 
metallireducensei_268'i). These two genes were added to the recon- 
struction by associating them to the appropriate reaction. In the case 
of the 23 genes where a discrepancy existed between the updated 
genome annotation and ModelSEED annotation, the updated 
annotation was taken as gold standard. No reactions were added or 
removed from the reconstruction as a result of this analysis. 



Table 3. Phenomic and modeling data from growth screens of G. metallireducens GS-15 wild type and dcuB with acetate as an 
electron donor. 



Donor Uptake 



Strain 


Acceptor/Growth Mode 
(Donor Acetate) 


Growth Rate 
(hr-1) 


Rate (mmol 
gDW-1 hr-1) 


Acceptor Uptake Rate 
(mmol gDW-1 hr-1) 


Acceptor/Donor Ratio 


GS-15 dcuB 


Fumarate/Batch 


0.114±0.005 


6.81 ±0.10 


14.0±0.2 


2.06±0.01 


GS-15 dcuB 


Fumarate/Cherrtostat* 


0.05 


2.79±0.23 


N.D. 


N.D. 


GS-15 


Nitrate/Batch^ (byproduct ammonia) 


0.1 11 ±0.022 


5.20±2.08 


4.67±0.64 (no3) ~ 
3.98±0.38 (nh4) 


0.90±0.42 


GS-15 


Fe(lll)/Batch 


N.D. 


N.D. 


N.D. 


6.96±0.43 


GS-15 


Fe{lll)/Chemostat» 


0.05 


7.86 ±1.02 


58.29±9.08 


7.41 ±0.64 


Simulation GS-15 

dcuB 


Max Growth Rate Fumarate 


0.116 


6.22 


14.2 (lim) 


2.28 


Simulation GS-15 

dcuB 


Max Growth Rate Fumarate 


0.054 


3.02 (Mm) 


7.06 


2.34 


Simulation GS-15 


Max Growth Rate Nitrate 


0.118 


7.28 (lim) 


4.35 -3.61 (lim) 


0.60 


Simulation GS-15 


Max Growth Rate Fe(lll) 


0.047 


8.88 (lim) 


62.31 


7.02 



* Donor-limited chemostat; 

N.D., Not Determined; lim, limiting rate in the simulation; 

'^calculated error using a 90% confidence interval from Lovley and Phillips, 1988 [27]. 
doirl 0.1 371/journal.pcbl.l003575.t003 
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Generation of the Biomass Objective Function (BOF) 

The BOFs were formulated using a previous template [36] and 
are included in Dataset S3. The biomass content previously 
determined for the close species Geobacter sulfuneducens was used to 
determine the breakdown of macromolecules [12] except that for 
total carbohydrate as the distribution in the murein, lipopolysac- 
charide, and cytosolic fractions was not indicated. Further, the 
genome annotation [14] was used for the breakdown of 
chromosome bases, a study on lipid and lipopolysaccharide chain 
length was used for the breakdown of acyl chain length [37], and 
the remaining content was approximated using the full profile 
presented for the gram-negative bacterium E. coli [36]. It should be 
noted that prediction of growth rate and unmeasured uptake rates 
are relatively insensitive realistic variations in biomass macromo- 
lecular weight fractions [36]. The BOF components included in 
the core BOF were extrapolated from the core BOF formulated 
for E. coli as Geobacter sp. have the same gram-negative cell 
morphology. 

Flux balance analysis sinnulations 

The COBRA Toolbox 2.0 [38] and the SimPheny framework 
(Genomatica, Inc., San Diego, CA) were used for simulations. 
Constraints used to simulate growth in ferric citrate medium are 
presented in Text SI. For the evaluation of the maintenance 
energies, the best-fit values of acetate and fumarate were set to the 
exact values. 

In vivo growth screens 

Growth screens were performed in triplicate using cultures in 
125 mL serum botdes under anoxic conditions using Fe(III) or 
fumarate (with the dcuB mutant studies) as an electron acceptor. 
The composition of the fumarate medium and ferric citrate 
medium were the same as previously described [39] and [40], 
respectively. The concentrations of the electron donors utilized in 
the experiments were 15 mM acetate, 20 mM ethanol, and 
10 mM butanol. For the growth screens, the cells were passed in 
the media in which they were being tested for at least three 
passages before the growth screen was performed. For the dcuB 
mutant strain, the correlation between OD and biomass that was 
used was 0.4561 gDW L-1 OD600-1. This value was determined 
by growing the dcuB cells in freshwater medium, taking OD 
measurements at various time points, and weighing the dried 
biomass at various points of the growth curve. The cells were dried 
overnight in an oven and weighed on filter paper, with a 
correction for the amount of weight lost by a filter paper that did 
not contain cells. Analytes were quantified by HPLC using an 
Aminex 87-H ion exchange column at 65° C. The mobile phase 
was 5 mM H2SO4 at an isocratic flow of 0.5 mL/minute. Sample 
injection volume was 10 \xL. Products were identified by retention 
time using ultraviolet detection at 210 nm and refractive index 
detection at 30° C internal temperature and 45° C external 
temperature and quantified by relating peak area to those of 
standards. Fe(II) was measured using the ferrozine assay as 
described [40]. Growth rates were calculated by determining the 
exponential growth phase region from a series of samples taken for 
each growth screen (typically, 7-10 samples were taken over the 
entire screen). For this corresponding exponential growth phase 
region, the ratio of an analyte to the gDW of the sample was 
determined using a linear fit obtained through the least-srjuares 
method ('regress' function in MATLAB). This value (mmol gDW- 
1) was then multiplied by the growth rate to get a corresponding 
uptake or production rate. The ratios were subsequendy 
multiplied by the growth rate to get the uptake or production 
rates. Averages and standard deviations were reported and 



calculated from three biological replicates for each experiment. 
For the G. metallireducens GS-15 dcuB stain, the fiimarate and 
succinate rates were independently calculated and then averaged 
for each replicate as there is a 1:1 ratio for the dcuB exchange and 
experimental accuracy for each analyte differed slightly. The 
fumarate uptake rate was calculated using the net fumarate and 
malate concentrations measured [41]. Malate was observed when 
growing with ethanol and butanol as electron donors with the dcuB 
strain. For the acceptor rates in these two conditions, the succinate 
production rate was used as the rate of the donor to minimize 
compounding measurement error (although the effective fumarate 
uptake rate was very similar, <27% difference in any individual 
replicate). 

Chemostat cultures of G. metallireducens GS-15 wild-type and 
dcuB cultures were grown at 30°C in anaerobic continuous culture 
vessels as previously dc-scribixl [42]. Tlu; GS-15 wild-t\pe strain 
was grown under donor limiting conditions with 5 mM acetate 
and 55 mM ferric citrate, the dcuB strain was grown under acetate 
limiting conditions with 5 mM Acetate and 27.5 mM fumarate. A 
protein conversion value of protein content (ug/ mL) = 0.4344*dry- 
cell-weight (ug/mL) was used to calculate rates. 

To determine if G. metallireducens was capable of growth with 
formate, cells were adapted to 20 mM formate in a ferric-citrate 
medium [40] and not provided with any other electron donor. 
After the third transfer with formate as the sole electron donor, 
samples were withdrawn anaerobically and substrate consumption 
was monitored. Formate was detected via high performance liquid 
chromatography. Reduction of Fe(III) to Fe(II) was monitored 
with the ferrozine assay. Experiments were run with triplicate 
cultures. The number of cells was monitored over time until the 
cultures reached a plateau of Fe(II) production. To fix cells, 
900 |a.L culture was withdrawn anaerobically and mixed with 
100 ixL glutaraldehyde (25%). Fixed cells were stored at -20°C 
until processed. All samples were defrosted, placed on filters, 
filtered and washed with sterile water prior to staining with acridin 
orange for 3 minutes. Filters were washed and dried. Cell counts 
were carried out for three biological samples by manually counting 
at least 5 microscopic fields (Area 5826 (tm^ per field) per sample, 
with the help of a cell counting application in ImageJ (http:// 
rsbweb.nih.gov/ij/). Control cultures consisting of (i) cells in ferric- 
citrate medium with no formate, and (ii) ferric-citrate medium 
with formate and no cells (i.e., cell-free) were also performed. 

Transcriptomic analysis 

For the microarray experiments performed, G. metallireducens was 
grown anaerobically with Fe(III) citrate (55 mM) as the electron 
acceptor and one of the following electron donors: acetate 
(10 mM), benzoate (1 mM), phenol (0.5 mM), or toluene 
(0.5 mM). Cells were grown in 1 L botdes and harvested during 
mid-exponential phase by centrifugation. The cell pellet was 
immediately frozen in liquid nitrogen and stored at —80 °C. RNA 
was isolated from triplicate cultures grown on each electron donor 
with a modification of the pre\ iousiy described method [43]. 
Briefly, cell pellets were resuspended in HG extraction buffer [44] 
pre-heated at 65 °C. The suspension was incubated for 10 minutes 
at 65 °C to lyse the cells. Nucleic acids were isolated with a phenol- 
chloroform extraction followed by ethanol precipitation. The 
pellet was washed twice with 70% ethanol, dried, and resuspended 
in sterile diethylpyrocarbonate-treated water. RNA was then 
purified with the RNA Clean-Up kit (Qiagen, Valencia, CA, USA) 
and treated with DNA-free DNASE (Ambion, Woodward, TX, 
USA). The RNA samples were tested for genomic DNA 
contamination by PGR amplification of the 16S rRNA gene. 
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cDNA was generated with the TransPlex Whole Transcriptome 
AmpUfication Kit (Sigma). 

Whole-genome microarray hybridizations were carried out by 
Roche NimbleGen, Inc. (Madison, WI, USA). Triplicate biolog- 
ical and technical replicates were conducted for all microarray 
analyses. Cy3-labeled cDNA was hybridized to oligonucleotide 
microarrays based on the G. metallireducens genome and resident 
plasmid sequences (accession number NC007515 and NC007517 
at GenBank). All the microarray data has been deposited with 
NCBI GEO under accession number GSE33794. 

For each metabolic shift (benzoate vs acetate, toluene vs acetate, 
phenol), the fold-changes in expression level and p-value (t-test) 
were computed using ArrayStar 4.02 (DNASTAR, Madison, WI, 
USA). This was used in conjunction with the iAF987 metaboUc 
model to predict the metabolic adjustment using the MADE 
algorithm [21]. Additionally, the constraints for the L\F987 
metabolic model were set to simulate growth on the respective 
substrates. MADE analysis was implemented using the TIGER 
toolbox [45]. The output of this analysis consisted of the genes 
predicted by MADE algorithm to significantly change in 
expression during the concerned metabolic shift. 

Supporting Information 

Dataset SI An Excel file containing the updated genome 

annotation for G. metallireducens. 

PCLSX) 

Dataset S2 An Excel file containing the content of the iAF987. 
pCLSX) 

Dataset S3 The biomass objective function formulated for G. 

metallireducens and a table corresponding to changes which were 
made to the network based on functional testing using the biomass 
objective function. 
(XLSX) 

Dataset S4 SBML files of the iAF987 reconstruction. 
(ZIP) 

Dataset S5 Transcriptomic data for G. metallireducens on acetate, 
benzoate, toluene, and phenol as electron donors and Fe(III) as the 
electron acceptor. 
(XLSX) 

Figure SI A Phenotypic Phase Plane (PhPP) analysis to 
approximate the energetic cost of external electron transfer. 
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